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Abstract 

Starting from the Flohlich electron-phonon model, which has a history 
of 50 years, a new unitary transformation is proposed to implement the 
perturbation treatment. Our main results are: (l)The phonon-induced 
interaction shows a crossover from the BCS-like potential when the phonon 
frequency ujp is much smaller than the Fermi energy Ep to that of the 
small polarons when uOp/ Ep > 1. (2)The jump of momentum distribution of 
electron number nk at the Fermi surface goes to zero when the dimensionless 
coupling constant A increases to the critical value Ac > 1, which means a 
possible broken down of the Fermi-liquid description. 



PACS numbers: 71.38.+i; 74.20.-z 
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It is a well-known fact that within the Migdal-Eliashberg (ME) description |1|, ^ 
of electrons and phonons coupled by the linear electron-phonon interaction there is no 
instability (where the Fermi-liquid description may break down) at any value of the 
dimensionless coupling constant A. People tried to study the polaronic collapse of the 
electron band starting from the Lang-Firsov (LF) transformation]^ followed by the 
small polaron approximation |P, 0. Note that the LF transformation together with the 
small polaron approximation cannot lead to those results which can be obtained via 
the ME approach. For example, Bardeen, Cooper, and Schrieffer[|] (BCS) proposed 
a square-well potential for electrons in momentum space, that is, only those electrons 
within a layer of width ujp (the characteristic phonon frequency) near the Fermi surface 
can attract with each other. This potential is localized in momentum space, so it is 
extended in real space and leads to large coherence length of Cooper pairs. But the 
LF transformation results in an attractive potential for all electrons in Fermi sea0, 0, 
which is a localized potential in real space. As far as we know, there is no theories 
which can describe successfully the crossover between the two pictures. In this letter we 
propose a new approach which can (l)describe the crossover and (2)lead to a possible 
broken down of the Fermi-liquid description. 

The Frohlich Hamiltonian of electron-phonon coupling system is|0 

H = 5I(ek - fio)dl^dk,a + ^q^i^q + ^ H H ^q4+q,<7C^k,a(&lq + &q)- (1) 
k,cr q VJV q 

is the number of sites, Ck is the bare band function. 6q and (ik,o- are usual notations 
for phonon or electron operators, respectively, /iq is the chemical potential, cUq the 
phonon frequency and Qq^ the electron-phonon coupling. We set h = 1 and fc^ = 1. 

Frohlich used a unitary transformation to treat H^, H' = exp{S)H exp{—S). The 
transformation can proceed order by order, 

H' = H, + H^ + \S, Ho] + [S, Hr] + ^[5, [S, H^]] + 0(4), 

where H = Hq+Hi and Hq contains the first two terms and Hi the last one. Frohlich let 
Hi + [S, Hq] = to get the generator S. The transformation leads to a phonon-induced 
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interaction with the potential 

2 

VF{k + ci,k) = - ^-^^ -, (2) 

|ek+q - CkT - t^q 

which can be attarctive or repulsive with a singularity at the energy shell: |ek+q — ek| = 
cUq. Frohlich noted that, although the transformation can eliminate the first order terms 
completely, because of the singularity one has to set a constraint | |ek+q — ek| — cjql > e 
(e > is a constant) in the transformation 0. Because of the constraint the elimination 
of the first-order terms is not conplete. 

The BCS theory simplifies the Frohlich potential as|^ 

T/ n , IN i -g'^/^p, for|ek+q - Ckl < 
V^C5(k + q,k) = | otherwL. 

We propose to improve the Frohlich transformation by a new generator 

^ = ^ E E -4+q,.rfk,.(&^-q - 6q)5(k + q, k) , (4) 

where (5(k + q, k) is a function of the energies of incoming and outgoing electrons in 
the electron-phonon scattering process, 

5{k + q, k) = 5(6k+q, 6k; cu^) = (^1 + . (5) 

The reason of choosing this generator will become clear lator. After transformation 
the first order terms in H' are 

Hn = H, + [S, i^o] = ^ E E -—JT -A^n,^^^^'^ 

q k,<7 ^q + l^k-q - Ckl 
X {[|ek-q - Ckl - (Ck-q - ek)]&-q + [jck-q " Ck | + (Ck-q - ek)]&q} • (6) 

The second order terms in H' are: 
Hn = [S,H,] + ]^[S, [S,H,]] 

= ^ E E ^(^'-q^-q + &q&l)^'(k + k) ^^^^^^ (4+q..4+q,. ' 4..4,. 

k.q q q 
E E -^(k + q, k)[2 - 5(k' - q, k')]4+q,.4,.4'-q,.'4',.', (7) 

k,k',qa,(T' ^q 



where the non-diagonal phonon terms are neglected. H12 contains a phonon-induced 
interaction with potential 

V{k + q, k) = -^5{k + q, k) [2 - S{k, k + q)] . (8) 

Fig.l shows y(k', k) as functions of ek' — Ck for different ratios ujp/ Ep where Ep is the 
Fermi energy. For comparison we also show VBcs(k', k), which is a narrow square- well 
in the middle, and the potential for small polarons which can be obtained via the LF 
transformation: Vlf — —g'^/^p for all electrons. 

The purpose of our transformation is to find a better way to divide the Hamiltonian 
into the unperturbed part and the perturbation. After transformation H' = H'q + Hji + 
Hj2 + O(giq). The unperturbed part is 

^0 = E ^q^i^q + E(^k - /X0)4,aC^k,a " E ^(k) (4,1^^^,^ + d-^,idu,^) , (9) 
q k,(T k,iT 

i^k = 6k + ^ E ^ COtH^f-^±^6\k + q, k) 

q ^q ^q 

E ^^(k + q, k)[2 - 6{k + q, k)]%^ tanh |±^, (10) 

iV q ^k+q 

A(k) = ^ E ^'^(k + q, k) [2 - 5(k + q, k)]^^±^ tanh |±a, (11) 

I\ q ?k+q -^-t 

Ck = V(^k-//o)^ + A2(k), (12) 

which can be solved exactly. The perturbation is Hji + Hj2 where Hj2 — Hj2 — 
{Hq — Ho). £'k and A(k) have been determined by the condition that the lowest order 
contribution of Hj2 to the self-energy is zero. 

A renormalized chemical potential /j, and a band renormalization factor p(ek) can 
be introduced, E]^ — fj,Q = p(ek)(ek — /x). The equation to determine /i is 



1 p(ek)(6k-/.) 1^ 

£k 2r' 



1 — n = — tanh 

AT ^-^ C. 
k 

where n is the electron number density. The equation to determine p(ek) will be given 
later. 
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The Green's function of is[p 



Go(k, uj) = {u; + E^- /io)/(cu^ - F^{k, u) = -A{k)/{u' - ^'(k)) 

The contribution of Hn to the self-energy (to the second order of gq) iscitemah 

,2 
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EG'o(k-q, 



UJ — ILOn 



Ck-q — Ckl — (Ck-q — Ck)]^ [l^k-q — ^kl + (^k-q " Ck)]^ 



iy(k,cj) = 0. 



(13) 
(14) 



The abnormal self-energy W{k, cu) = because it contains a factor |ek-q — Ckp — (ck-q — 
ek)^ in the q-summation. Generally, the normal self-energy S(k, cj) 7^ 0. But for the 
normal state (A(k) = 0) and when T = 0, we have 



Sn(ek = Ai,^^) = 0, 



(15) 



that is, the normal self-energy is zero at the Fermi surface. As the spectrum of ele- 
mentary excitations in the normal state is = p(ek-q)(ek-q — /i) + Sn(k, tu), the mass 



renormalization at Fermi surface is[|TT 

Tfl 3 

— = P(ek) + ^S„(k,^) 
m* ack 



d 

1 - — S„(k,u;) 



p(ek = /i). (16) 



Eqs.(14), (15), and (16) are main reasons for the choice of the functional form of 
(5(k + q, k) in (5). Now one can see clearly the purpose of our unitary transformation: 
The transformed H' is divided into the unperturbed part H^, which contains the main 
physics of the problem, and the perturbation Hji + Hj2, which is small as shown by 
Eqs.(14), (15), and (16). 

Gap equation (11) can be rewritten as an integral equation by introducing the 
Eliashberg function a^-F^, ^: 



' de' r.l^^.-F((^),^ + f'-^-:'^tanhi(!:i. (17) 
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N{e) is the density of states (DOS). The integration over e' is from the bottom (— -D) 
of the band to the top {D). The equation for p(ek) can be rewritten in the same way, 



rD POO Nff') O 

p(e) = 1 + / de' c/^]^^a'F(fi) coth( 
J-D Jo N(ii) 



fx 



(fi+|e'-e|)2 (n+\e'-i^\yj 



ek-/i) 



f de' rdQ-^a'F{n) , 

7-D 7o A^(/i) 1 (f^ 



e\)^ 



n + 2\e' 




(n + \e' - 





2T 



These two equations are very similar to the Ehashberg equations[0, ^ 

For calculating the physical quantities we must calculate the thermodynamical po- 
tential and the average of electron or phonon operators. The thermodynamical poten- 
tial is 

n = -^lnTrexp[-/3/7] = InTr exp[-/3if'] ^ InTr exp[-/3if;]. (19) 

fj fj jj 

The last "~" is because of Eqs.(14) and (15), that is, to the order 0{g^) the contribution 
of Hji + Hj2 is very small in the lowest temperature region. Hence, for the normal 

state 

2 

Q ^ -— ^ln{l + exp[-/5p(ek)(ek - /i)]} + phonon part. (20) 
pjv k 

The heat capacity can be calculated as C = —Td'^fl/dT^ = Co/p(ek = p), where Co is 
the heat capacity for free electrons and 

In calculating we assume a constant DOS and A = 2/ dQa'^F{Q)/Q. The enhancement 
factor at ujp/D is 1/(1 — A) which should be compared with the same factor in 
ME theory 1 + A[|. 

We have to take into account the effect of the unitary transformation when cal- 
culating the thermodynamical average of electron or phonon operators. The Green's 
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function for the original Hamiltonian H 



G(k, r) = -Tr {t, exp[-/3(i7 - r])]rfk,.(r)<,} 

= -Tr {T^exp[-f3{H' - Q)K^Ar)dt} , (22) 



where means r ordering[|rT|]. The transformation of a single fermion operator can 
proceed as 

dk,a = e^4,<xe-^ = 4,<x + [S, rfk,a] + ^[5, [S, rfk,.]] + 0(4). 

/3 



G(k,uj) = I rfrG(k,r)e' 

JO 



/o 

= ^(k, ^) - ^ E k - q)G'o(k, u) 

E k - q)^ E — ^<^o(k - q, ^ - iw„), 

where G(k, u) is the Green's function of H'. Starting from the Green's function one 
can calculate various physical quantities. As an example, we calculate the number of 



electrons in a momentum state k for the normal state when T = 0[|11 
1 fO 



rik = / dujlmG{k, u + iO"*") 

TV J-oo 

= f 1 - ^ E 4^'(k> k - q)) ^(/x - ek) + ^ E 4'^'(k, k - q)e(/i - ek-q). 
\ q q / q q 

Fig. 2 shows rik as functions of ek around the Fermi surface. For smaller coupling 
(dashed line) or the larger frequency (dash-dotted line), there is a finite jump of nk 
at Fermi surface. But for smaller frequency (solid line) and lambda ~ 1, the Fermi 
surface is smeared by the electron-phonon coupling. The jump of nk at Fermi surface 
ek = /i is p(ek = yu) (Eq.(25)) and it predicts an instability of the Fermi liquid state at 

A.=(i- y . (23) 

Ac — > 1 when cUp — > 0. For comparison, ME theory predicts a jump 1/(1 + A) and there 
is no instability. 
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At the end, as the transformation is truncated after the second order of g^, we 
justify the cutoff by showing the small expansion parameter. Roughly speaking, in 
three-dimension it is Xup/Ep when Up/Ep < 1 but XEp/up when ujp/Ep > Icitezh. 
Note that they are the same as that of ME theory or LF transformation. 
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Figure Captions 



Fig.l The phonon- induced interaction. 

Fig. 2 The phonon- induced renormahzation of the momentum distribution on 
and near the Fermi surface. 
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